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Summary 

A weak-wave analysis of shock interaction with 
a slipstream is presented. The theory is compared 
with the acoustic theory and to the exact nonlinear 
analysis. Sample calculations indicate that the weak- 
wave theory yields a good approximation to the exact 
solution when the shock waves are sufficiently weak 
that the associated entropy increase is negligible. A 
qualitative discussion of the case of counterflowing 
streams is also included. 

Introduction 


R reflection coefficient, defined by 

equation (35) 

T transmission coefficient, defined by 

equation (36) 

z parameter defined by equation (10) 

f) = VM 2 - I 

'y ratio of specific heats 

A change in quantity due to shock or 

expansion wave 

8 small variation in quantity 


When a shock is incident on an interface of rela- 
tive motion, a wave is usually transmitted through 
the interface and another wave is reflected. The 
strengths of these two waves are determined by two 
boundary conditions: pressure and flow directions 
are equal on the two sides of the interface both ahead 
of, and behind, the interaction point. The corre- 
sponding problem for sound waves has been solved 
by both Ribner (ref. 1) and Miles (ref. 2). These 
solutions corrected earlier erroneous analyses. For 
shock waves, the changes in flow direction are re- 
lated to wave strengths in a highly nonlinear way; 
consequently, solutions are obtained by iteration (as, 
for example, in ref. 3). 

The weak- wave analysis, which is given herein, 
provides an intermediate solution valid for weak 
shocks. It yields analytic expressions which are more 
accurate for shocks than the acoustic approximation. 
Furthermore, the derivation of the results provides 
some physical insight into the relationship of the 
acoustic theory to the full shock wave problem. 

The material is presented in the following order. 
First, the problem is described and the boundary 
conditions at the interface are given. Then, the 
results of references 1 and 2 for the acoustic case are 
summarized. The acoustic results are taken as a first- 
order approximation for weak shocks, and the weak- 
wave theory is developed as a correction to this first- 
order approximation. One way of solving the exact 
equations is then described and the results of the 
three theories are compared. Finally, a qualitative 
discussion of the case of counterflowing streams is 
given. 

Symbols 

a coefficient defined by equation (17b) 


£ 

0 

v 

Subscripts: 

i 

r 

T 

t 


o 

1,23,4,5 


= ^ 
p 

shock angle 

flow deflection angle 

incident 

reflected 

total 

transmitted 
acoustic solution 

region 1,2, 3, 4, 5, respectively (fig. 1) 


Analysis 

Basic Considerations and Boundary Conditions 

If the incident shock is curved, its radius of cur- 
vature is assumed to be large relative to the thick- 
nesses of the shock waves and of the slipstream. Con- 
sequently, the interaction may be treated as a local 
phenomenon (assuming negligible shock wave thick- 
ness). (See, for example, the interesting schlierens in 
ref. 4.) 

The flow configuration is diagrammed in figure 1. 
The undisturbed parallel flows in regions 1 and 5 
proceed from left to right. A shock i in region 1 is 
incident on the slipstream, giving rise to a reflected 
compression or expansion wave r and a transmitted 
shock t. 

The pressures in regions 1 and 5 are equal, as are 
those in regions 3 and 4. Consequently, the pressure 
increments through the waves satisfy the relation 


A p % + Ap r = (p2 — pi ) + (p3 — P2) — P 4 - Ps — Apt ( 1 ) 


b coefficient defined by equation (17c) 

M Mach number 

p pressure 


Similarly, the changes in flow direction must 
match on the two sides of the slipstream. Notice, 
however, that the incident shock deflects the flow 
downward, whereas a reflected shock tends to deflect 


the flow upward. Therefore, this latter deflection ap- 
pears with a negative sign. 


Vi~v r = v t (2) 

Thus, a reflection expansion produces a negative u r . 

Equations (1) and (2) are exact within the 
assumptions of the theory. Approximations en- 
ter in converting the pressure increments to shock 
strengths and in relating the shock strengths to flow 
deflections. 


where 

75 ^/ 3 , 


The solution of this set of equations is 


R = 
T = 


1 - z 
1+2 

2 

1+2 


( 10 ) 


( 11 ) 

( 12 ) 


Acoustic Theory 

As was mentioned in the Introduction, the acous- 
tic approximation has been treated independently by 
Ribner (ref. 1) and Miles (ref. 2). Ribner’s treat- 
ment is especially interesting in the present context, 
because he treats the sound waves as Mach waves. 
However, in the usual acoustic analysis the pressures 
P2, P 3 > and p 4 represent sinusoidal perturbations on 
the undisturbed pressure p x (= p 5 ), and equation (2) 
represents a relation between the slopes in a sinu- 
soidal ripple in the slipstream. 

Equation (1) is written in the form 


Ap, A pr 

Pi Pi 


Ap t 

P 5 


In the acoustic approximation, 

P Ap 

7M 2 p 


v — 


( 3 ) 


( 4 ) 


Consequently, in this approximation, equation (2) 
can be written 


Pi f &Pi _ Ap, 

71 Mj 2 V Pi 


P5 Apt 


Pi ) 15 M? p 5 


( 5 ) 


If the reflection coefficient is defined by 


This acoustic solution represents a first-order so- 
lution for weak shocks. Since it is an acoustic ap- 
proximation, it is independent of wave strength. It 
also provides other useful information relating to the 
shock wave solutions. For example, if 2 = 1, then 
R = 0 so that there is no reflection, and the entire 
wave is transmitted. One such case occurs when the 
two streams are at the same temperature, there is no 
relative motion, and 7 j =75. But it can also occur 
if there exists a small relative motion provided that 
there is a compensating difference in the 7’s. . 

Equation (11) indicates that when z < 1 the 
reflected wave is a compression, and equation (12) 
indicates that the transmitted wave is stronger than 
the incident wave. Conversely, for 2 > 1, the reflected 
wave is an expansion and the transmitted wave is 
weaker than the incident wave. These qualitative 
results should be applicable within limits for shock 
waves. They are useful in selecting the appropriate 
set of equations to be solved for full-shock solutions 
and also in setting the limits for the intervals over 
which solutions are sought. 

Weak-Shock Theory 

To treat waves of finite strength, equation (1) 
must be expressed in terms of wave strengths. This 
is accomplished by writing the second term in equa- 
tion ( 3 ) (which is an exact equation) in the form 


R = Ap r /p 2 _ A pr/pi _ Apr 
A Pi/pi A pi/pi ~ Ap, 

and the transmission coefficient by 


( 6 ) 


^Pr _ A p^ P2 _ A p r / P 2 — pi \ 

Pi P2 Pi P2 V Pl ) 

Equation ( 3 ) then becomes 


— £ r(l + £,') 

( 13 ) 


T _ &Pt/pj> _ A pt/pi _ Apt 
A Pi/Pi Api/pi A pi 


( 7 ) 


£ i + (1 + £ i) £ r = £t ( 14 ) 


equations ( 3 ) and ( 5 ) become, respectively, 

1 + R = T (8) 

1 - R = zT (9) 


It is assumed that either the incident shock strength 
or the incident shock angle is prescribed, since each 
of these quantities can be directly prescribed in terms 
of the other. (See ref. 5 , formula ( 128 ); see also 
subsequent eq. ( 26 ).) 
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The weak- wave solution can be treated as a cor- 
rection to the acoustic solution. Thus, denoting the 
results of the acoustic solution (eq. (4)) by the sub- 
script o, equation (2) becomes 


terms higher than second-order yields 

+ b\v\ + a 2 V r ,o + a 2 bv r + a l a 2 v i v r,o 

+ b 2 — a 5 "t,o + 05 + h^t,o ( 20 ) 


— (Vr,o d" bl^r) — ^t,o bui ( 16 ) 


which, by subtracting the acoustic relation, yields 
simply 

— 8i/ r = But (16) 

The second-order relation between the flow de- 
flection angle and the wave strength is (see ref. 5, 
formulas (151) and (174)) 


where 


e = av + bu 2 

(17a) 

' yM 2 

P 

(17b) 

b=^[^+l)M 4 -4d 2 } 

(17c) 


The values of 7 and (3 in these expressions are 
those immediately ahead of the wave. For an ex- 
pansion, the deflection angle v is negative and, con- 
sequently, the first-order term in equation (17a) be- 
comes negative. The second-order term is identical 
for compression and expansion waves. 

Since 71, 75, and M5 are prescribed, b\, 
05, and 65 can be computed directly. In order to 
compute a 2 and 62, the value of M 2 is required. It 
could be estimated by the weak-wave approximation 
(ref. 6, p. 292), but since the exact expression (ref. 5, 
formula (157)) is only slightly more complicated, it 
was used in the present study. 

Substituting equations (17) into equation (14) 
yields 

a\V{ + b\uf + (l + °1 v i + ( a 2 v r + ) 

= a^t J t-b^t (18) 

If terms higher than second-order are discarded, 
equation (18) becomes 

Q L 2 

a\v^ + b \u i + a 2 v r + 01 a 2 ViV r + o 2 v r 

= 051^ + 65^ (19) 

Now when v is written as a correction to the 
acoustic solution, as in equation (15), the correction 
8v results from including terms of order v 2 . Thus, 
substituting v — v Q 4~ bv for the reflected and trans- 
mitted waves in equation (19) and again eliminating 


Subtrat ting the first-order accurate relation 

a\ Vi + a 2 ^ r ,o = ( 21 ) 

yields for the second-order terms 

2 2 l 2 

■12b Vr - a$8vt — b^v to - b 2 v r , 0 - b \v { 

- aia 2 ViV rf0 (22) 


The interpretation of the terms in this equation is 
as follows. The coefficients a 2 and <25 of the left- 
hand terms indicate the proportional influence of the 
second order corrections on the reflected and trans- 
mitted wave strengths, respectively. The first three 
terms on the right result from including the nonlinear 
second order term in the relation between flow angle 
and wave strengths. The last term represents the in- 
fluence of the pressure increase through the incident 
wave on the strength of the reflected wave. This lat- 
ter term is negative for a shock and positive for an 
expansion. 

Now substituting — 8u r for 8vt from equation (16) 
into equation (22) yields the solution 


8v r 


h^t,p ~ b 2^r y o ~ hVj ~ a l a 2 ViVr,o 
02 + 05 


(23) 


The actual flow deflection through the reflected wave 

v r = Vr } o + bv r (24) 

The flow deflection through the transmitted wave can 
now be obtained from equation (2). The strength of 
the reflected wave is, from equation (17a), 


Sr — CL 2 V r + b 2 V 2 


(25) 


and St is then obtained from equation (14). 

Exact Shock Relations 

To solve the exact equations for the flow quan- 
tities n regions 3 and 4 requires a numerical pro- 
cedure. The specific formulas used in the procedure 
depend on the nature of the reflected wave (compres- 
sion or expansion) and also on the flow parameters 
that are specified. Any fundamental quantity relat- 
ing to the reflected or transmitted waves or to the 
flow in regions 3 or 4 may be chosen as the unknown 
parameter to be solved for in the numerical proce- 
dure. For the following calculations, e r was chosen 
for the unknown parameter if the reflected wave was 
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a compression, but M3 was found to be a more con- 
venient parameter to solve for if the reflected wave 
was an expansion. 

Reflected compression. For the former case (re- 
flected compression) the following method was used. 
The quantities Mj, M 5 , qq, 75, Pl) p 5 , and either 
£i or 0 t are assumed to be given. If 9 t rather than 
£i is specified, then e t is computed from the formula 
(ref. 5, formula (128)) 


angles now, since they are computed at each step 
in the numerical procedure and written over at the 
subsequent step. Thus, the final values stored in 
these locations when the procedure has converged are 
the correct values. 

Reflected expansion. When the reflected wave 
is an expansion, it is convenient to choose for the 
unknown parameter the Mach number in region 3 
M3; then the relation (ref. 5, formula (44)) 


27 M 2 sin 2 9 — (7 — 1) 
7+1 


(26) 


Note that here e corresponds to £ — 1 in reference 5. 
The Mach number in region 2 is obtained from the 
relation (ref. 5, formula (157)) 


M 2 = M 1 2 [(7+l) £ + 27l-2e(e + 2) 

2 (l + e)[( 7 - l)e + 2 7 ] 


(27) 


It is assumed that the value of 7 does not change 
through the shock (although different values of 7 may 
exist on the two sides of the slipstream). The flow 
deflection v x due to the incident shock is determined 
by (ref. 5, formula (160)) 


(_i_ ) J MMl - 1) - i£t Ilf 

V7A/ 2 - e ) y (7 + l)e + 27 

^ 1 . ( 2g ) 

Thus, the required quantities in region 2 are com- 
puted directly without iteration. Equation (28) also 
yields a relation 


v r = V r (M2,€ r ) (29) 


P_ 

PT 


1 + 


7 - 1 


M 4 




(32) 


is applied as follows. First, using the known values of 
P2 and M 2 (eq. 27)), equation (32) is solved for pp 2- 
Then, since pp^ — pp,2 f° r a reflected expansion, 
equation (32) determines p 3 as a function of M3. The 
flow deflections, in expanding from M = 1.0 to M2 
and from M = 1.0 to M 3 , are calculated from the 
relation (ref. 5, formula (173b)) 


v — 



tan 1 


— cos 


1 

M 



(M 2 - 1) 


(33) 


The difference of these two calculations is v r (M 3 ), 
Now p 3 (M 3 ) is calculated by equation (32) and, since 
P2 is known, e r can be calculated. Substituting this 
value into equation (14) yields 67 (M3). The flow 
deflection v t associated with the transmitted shock is 
then calculated as a function of M 3 by equation (28). 
Finally, the value of M 3 is determined numerically as 
the zero of the function 


The transmitted wave strength et is now ex- 
pressed in terms of e r by equation (14). Equation 
(28) is then used to determine Vt = vt{M§,et). But 
with Et replaced by €t(e r ) and with 7 = 75, this re- 
lation becomes a function of e r 


g(M 3 ) = v { - v r (M 3 ) - v t (M 3 ) (34) 

The remaining flow variables can now be calculated 
by substituting this value of M 3 back into the previ- 
ous relations. 


V t = V t {M^6r) (30) 

Finally, substituting the functions in equations (29) 
and (30) into equation (2) yields a single equation 
for the unknown parameter e r , which is determined 
numerically as a zero of the function 

fi £ r) = v t - v r {e r ) - v t {e r ) (31) 

From the value of e r thereby obtained, e t is found 
from equation (14), and the reflected and transmitted 
shock angles are found by solving equation (26) for 
6. It is not necessary to determine the flow deflection 


Computed Examples 

Several sample calculations were performed for 
the purpose of comparing the acoustic, weak-wave, 
and exact solutions. The results of the calculations 
are displayed in figures 2 and 3. For the weak-wave 
and exact solutions, the reflection and transmission 
coefficients are defined by 


to 

III 

2} 

(35) 

rp 

~ Si 

(36) 
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Equations ( 6 ) and (7) for the acoustic case are con- 
sistent with these definitions since the acoustic pres- 
sure amplitudes are all referred to the undisturbed 
pressure pi (= P 5 ), and dividing both numerator and 
denominator by this quantity does not change the 
ratio. 

The results demonstrate that, for all the exam- 
ples, the weak-shock solution is an improvement over 
the acoustic approximation, and in most cases it is 
remarkably accurate. In all cases for which the re- 
flected wave is an expansion or a very weak shock, the 
weak-wave results are significantly better than those 
for a reflected shock of strength comparable with that 
of the incident shock. (Compare figs. 2(h) — (k).) 
This indicates that, although the entropy increase 
is third order in shock strength, the effect of the en- 
tropy variation is noticeable. The approximate theo- 
ries also yield poor results when the local Mach num- 
ber falls into the highly nonlinear region near M = 1. 
(See fig. 2 (a), where M 4 — > 1.) 


inconsistent, since large amplitude acoustic waves 
propagate nonlinearly and develop into shock waves, 
which are governed by a different set of boundary 
equations. For these equations one should assume 
that the free-stream conditions 71 , M\ , 74 , and M 4 
(fig. 4 (b)) are prescribed, as well as the strength £\ 
of the incident shock. However, at the interface, up- 
stream influences always occur since each stream is 
flowing upstream relative to the other. Thus, a flow 
angularity can develop ahead of both incident and re- 
flected shocks, and consequently, the problem is not 
well-posed. 

However, if one simply proceeds formally, the 
pressure relation at the interface becomes 

?2 ~ Pi + P3 ~ P2 P2 = / PS ~ P4 \ P± ^ 37a j 

Pl P2 PI V P4 / P5 

or 

£{ + (1+ ei)£ r = - - (37b) 

1 4 et 


Counterflowing Streams 

References 1 and 2 observe that strong amplifica- 
tions of the reflected and transmitted waves are possi- 
ble if the streams separated by the interface are flow- 
ing in opposite directions. It is difficult to imagine 
a practical steady-state situation for which two su- 
personic streams would be counterflowing. One can, 
however, conceive of several situations for which such 
a phenomenon might occur on a transient basis. One 
possibility would be a shock wave, emitted by a pass- 
ing airplane or from a blast, incident on a supersonic 
jet. Figure 4 depicts such a possibility. 

Suppose, for example, the jet emerges at M = 3 
relative to the ambient air, and an airplane is flying 
in the jet flow direction at M — 1.6 so that the 
shock that it generates is incident on the jet. Then, 
relative to a coordinate system moving with the point 
of interaction at the interface, there is a M = 1.6 
flow moving toward the left above the interface and 
(assuming the same sound speed in the jet as in the 
ambient air) a M = 1.4 flow moving toward the right 
below the interface. 

The acoustic analysis (refs. 1 and 2) of this prob- 
lem is straightforward. The Mach number of one 
of the flows is simply indicated to be negative and 
the parameter z in equation (9) is also determined 
to be negative. Then, the denominators in the ex- 
pressions for the reflection and transmission coeffi- 
cients (eqs. ( 11 ) and ( 12 )) can become very small 
and therefore yield large values for the reflected and 
transmitted wave amplitudes. However, it should 
be mentioned that this analysis rapidly becomes 


or, since 


P2 P3 
Pi P2 


= U + £i)U + £ r) — 


1 

1 4 St 


P4 

PS 


£{ 4(14 £i)z r - “(I4£i)(l4 S\ )£t (37c) 

Equation (37b) or (37c) may be compared with the 
simpler equation (14) for the coflowing case. 

A v/eak-wave analysis of this problem would 
hardly be justified inasmuch as the problem, as 
noted, s not well-posed, it is largely academic, and 
the weak-wave theory loses accuracy for strong re- 
flected shocks. 


Concluding Remarks 

A weak-wave analysis of shock interaction with a 
slipstream has been presented. The theory was com- 
pared with the acoustic theory and the exact nonlin- 
ear analysis. Sample calculations indicated that the 
weak- wave theory represented a good approximation 
to the rxact solution when the shock waves are suf- 
ficiently weak that the associated entropy increase 
is negligible. A qualitative discussion of the case of 
counter flowing streams was also included. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
September 19, 1988 
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Figure 1. Diagram of shock-slipstream interaction. 














(g) Mi = M 5 = 2.0; 71 = 1.2; 72 — 1*4 (z — 0.8571) 


Exact - 
Acoustic - 


Weak wave — 








(k) Mi = 1.5; M 5 = 4.0; 7l = 75 = 1.4 (z = 0.487) 
Figure 2. Concluded. 













(b) System of figure 4(a) in coordinate system moving with interaction point. 

Figure 4. Transient phenomenon giving rise to counterflow situation. Sound speed in jet assumed to be 
equal to that in ambient air. 
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